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LOCAL  HEAT  TRANSFER  AND  RECOVERY  TEMPERATURES  ON  A  YAWED 
CYLINDER  AT  A  MACH  NUMBER  OF  4.15  AND  HIGH  REYNOLDS  NUMBERS  1 

By  Ivan  E.  Beckwith  and  James  J.  Gallagher 


SUMMARY 

Local  heat  transfer ,  equilibrium  temperatures , 
and  wall  static  pressures  have  been  measured  on  a 
circular  cylinder  at  yaw  angles  of  0°,  10°,  20° ,  Ifi° , 
and  60°.  The  Reynolds  number  range  of  the  tests 
was  from  1X10 6  to  4'XIO6  based  on  cylinder 
diameter . 

Increasing  the  yaw  angle  from  0°  to  increased 
the  stagnation-line  heat-transfer  coefficients  by  100 
to  180  percent .  A  further  increase  in  yaw  angle  to 
60 0  caused  a  40-percent  decrease  in  stagnation-line 
heat-transfer  coefficients. 

At  zero  yaw  angle  the  boundary  layer  over  the 
entire  front  half  of  the  cylinder  was  laminar  but  at 
yaw  angles  of  40°  and  60°  it  was  evidently  completely 
turbulent ,  including  the  stagnation  line ,  as  deter¬ 
mined  by  comparison  of  local  heat-transfer  coefficients 
with  theoretical  predictions.  The  level  of  heating 
rates  and  the  nature  of  the  chordwise  distribution 
of  heat  transfer  indicated  that  a  flow  mechanism 
different  from  the  conventional  transitional  boundary 
layer  may  have  existed  at  the  intermediate  yaw  angles 
of  10°  and  20°.  At  all  yaw  angles  the  peak  heat- 
transfer  coefficient  occurred  at  the  stagnation  line 
and  the  chordwise*"  distribution  of  heat-transfer 
coefficient  decreased  monotonically  from  this  peak. 

The  average  heat-transfer  coefficients  over  the 
front  half  of  the  cylinder  are  in  agreement  with 
previous  data  for  a  comparable  Reynolds  number 
range. 

The  theoretical  heat-transfer  distributions  for 
both  laminar  and  turbulent  boundary  layers  are 
calculated  directly  from  simple  quadrature  formulas 
derived  in  the  present  report. 


INTRODUCTION 

Design  studies  of  hypersonic  lifting  vehicles  have 
generally  indicated  that  aerodynamic  heating  may 
be  reduced  by  using  highly  swept  configurations 
with  blunted  leading  edges.  For  laminar  boundary 
layers  the  effect  of  sweep  angle  A  on  the  heat 
transfer  at  the  leading  edge  is  usually  taken  as  cos 
A  as  shown  by  the  data  of  Feller  (ref.  1)  who 
measured  the  average  heat  transfer  on  the  front 
half  of  a  swept  cylinder.  More  recent  data  (refs. 
2  and  3)  have  indicated  that  the  effect  of  sweep 
may  be  more  nearly  cos3/2A  which,  at  a  sweep 
angle  of  75°,  would  result  in  a  50-percent  reduction 
of  the  heat  transfer  predicted  by  the  cos  A  varia¬ 
tion.  The  data  and  theory  of  reference  4  also 
indicate  a  cos 3/2  A  variation  but  the  theories  of 
references  5  and  6  indicate  a  variation  somewhere 
between  cos  A  and  cos3/2A  for  large  stream  Mach 
numbers. 

The  data  of  reference  7,  in  contrast  to  the  in¬ 
vestigations  just  cited,  showed  large  increases  in 
average  heat  transfer  to  a  circular  leading  edge 
with  increasing  A  up  to  a  A  of  about  40°.  These 
increases  in  heat  transfer  were  probably  caused 
by  transition  to  turbulent  flow  which  apparently 
resulted  primarily  from  the  inherent  instability 
of  the  three-dimensional  boundary-lay'er  flow  on 
a  yawed  cylinder.  The  leading-edge  Reynolds 
numbers  of  reference  7  were  .considerably  larger 
than  the  values  in  references  1  to  4  and  were  also 
larger  than  typical  values  for  full-scale  leading 
edges  of  hypersonic  vehicles;  hence,  the  main 
application  of  the  high  Reynolds  number  tests 
will  probably  be  to  bodies  at  angle  of  attack. 


i  Supersedes  declassified  NASA  Memorandum  2-27-59L  by  Ivan  E.  Beckwith  and  James  J.  Gallagher,  1959. 
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In  reference  8  it  was  shown  that,  even  for  bodies  of 
moderate  fineness  ratio  at  angles  of  attack  of  about 
20°  or  more,  the  heat  transfer  in  a  turbulent 
boundary  layer  is  nearly  the  same  as  that  for  a 
turbulent  boundary  layer  on  a  comparable  yawed 
cylinder  of  large  aspect  ratio.  The  laminar  theory 
of  reference  9  predicts  essentially  the  same  result; 
that  is,  the  heat  transfer  on  a  slender  cone  at  an 
angle  of  attack  of  about  15°  is  already  the  same  as 
that  on  a  corresponding  infinite  cylinder  at  75° 
yaw  or  sweep. 

The  present  investigation  was  undertaken  in 
an  attempt  to  understand  better  and  to  explain 
the  results  of  reference  7  by  means  of  local 
measurements  of  heat  transfer  to  a  circular 
cylinder.  The  tests  were  conducted  at  a  stream 
Mach  number  of  4.15  and  a  stream  Reynolds 
number  range  of  1X106  to  4X108  based  on 
cylinder  diameter.  The  experimental  results  are 
compared  with  theories  for  local  heat  transfer 
in  both  laminar  and  turbulent  boundary  layers  on 
yawed  infinite  cylinders.  The  laminar  theory 
is  based  on  the  similar  solutions  of  reference  10  for 
arbitrary  pressure  gradient.  The  turbulent  theory 
makes  use  of  the  integral  momentum  and  energy 
equations,  simplified  by  means  of  Stewartson’s 
transformation  (ref.  11),  together  with  suitable 
assumptions  for  velocity  and  thermal  profiles, 
skin  friction,  and  Reynolds  analogy. 
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SYMBOLS 


constant  coefficient  in  skin-friction 
law  (see  eq.  (3)) 
stagnation  speed  of  sound 
constant  in  equation  (BIO) 
specific  heat  at  constant  pressure 
cylinder  diameter 
integral  thickness  parameter, 


integral  thickness  parameter, 


stagnation  enthalpy,  cpT-\ 


u2+r2 

2 


heat-transfer  coefficient, 
heat-transfer  parameter, 


Q.w 

Taw-Tw 

T-Tw 

Taw—Tw 


e'w 


thermal  conductivity 


x 


Z 

z 

a 


7 

A 

8 

$* 


V 


reference  length 
Mach  number 

reciprocal  of  exponent  in  skin-friction 
law  (see  eq.  (3)) 

Prandtl  number  (0.7  used  in  all  cal¬ 
culations  except  as  noted), 
pressure 

resultant  velocity  in  transformed 
coordinate  system,  ^U2+  V2 
local  heat-transfer  rate  per  unit  area 
heat  transfer  in  transformed  coordi¬ 
nate  system  defined  by  equation 
(AS)  ■  ■ 

Reynolds  number  based  on  diameter, 

PUrD 

M 

recovery  factor 
temperature 


chordwise  velocity  in  transformed 
coordinate  system  (see  eq.  (A4)) 
chordwise  velocity 
resultant  velocity 

spanwise  velocity  in  transformed 
coordinate  s3~stem  (see  eq.  (A4)) 
spanwise  velocity 

chordwise  coordin  ate  (transformed 
system;  see  eq.  (A4)) 
chordwise  coordinate  (normal  to 
cylinder  axis;  physical  boundary- 
layer  system) 

coordinate  normal  to  surface  (trans¬ 
formed  system;  see  eq.  (A4)) 
coordinate  normal  to  surface  (physi¬ 
cal  boundary-layer  system) 
angle  between  local  streamline  and 
cylinder  axis 

pressure  gradient  parameter  (see  eq. 

(Bll)) 

ratio  of  specific  heats  (1.40  used  in  all 
calculations) 

thermal  boundary-layer  thickness 
velocity  boundary-layer  thickness  in 
transformed  coordinate  system 
integral  thickness  parameter, 

similarity  variable 
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Subscripts : 

av 

aw 

e 

L 

N 

o  1 

V 

r 

R 

s 


T 

t 

w 

X 

Y 


angular  distance  around  cylinder 
from  stagnation  line  in  a  plane 
normal  to  cylinder  axis 


enthalpy  profile  parameter, 


H-Hw 

Hx-Hw 


integral  thickness  parameter, 


.  f  KH£)> 

integral  thickness  parameter, 
integral  thickness  parameter, 


f 


(1  ~6)dZ 


yaw  angle  (complement  of  acute 
angle  between  free-stream  flow 
direction  and  cylinder  axis) 
viscosity  coefficient 
coefficient  of  kinematic  viscosity,  ju/p 
density 
shear  at  wall 

shear  in  transformed  coordinate  sys¬ 
tem  defined  by  equation  (A5) 


average 

zero  heat  transfer  at  wall  over  entire 
surface 

local  equilibrium  wall  conditions 
(vanishing  local  heat  transfer) 
laminar 

component  normal  to  cylinder  axis 
free-stream  stagnation  conditions 
external  surface  of  heat-meter  plug 
reference  point  within  boundary 
layer 
resultant 

value  at  stagnation  line  (when  used 
on  a  gas  property  such  as  p,  T, 
etc.,  denotes  that  the  property  is 
evaluated  outside  the  boundary 
layer  at  stagnation  line) 
turbulent 

stagnation  conditions  downstream  of 
normal  shock 
wall  conditions 
chordwise 
spanwise 

free-stream  conditions  ahead  of  bow 
shock 


1  local  conditions  just  outside  bound¬ 

ary  layer 

A  prime  denotes  differentiation  with  respect  to  ??. 

APPARATUS 

TUNNEL 

The  tests  were  conducted  in  one  of  the  blow¬ 
down  tunnels  of  the  Gas  Dynamics  Branch  at  the 
Langley  Besearch  Center.  The  tunnel  is  a  con¬ 
ventional  two-dimensional  nozzle  with  a  test  sec¬ 
tion  12  inches  wide  and  13  inches  high.  The 
nominal  Mach  number  in  the  test  section  is 
4. 15  ±0.03.  The  stagnation  temperature  for  these 
tests  ranged  between  150°  F  andJ300°  F.  Further 
details  of  tunnel  operation  and  flow  calibration 
may  be  found  in  references  7  and  12. 

MODELS  AND  INSTRUMENTATION 

Two  different  models  were  used  in  this  investi¬ 
gation.  The  models  were  circular  cylinders  made 
from  stainless  steel  and  were  the  same  in  all  re¬ 
spects  except  for  the  type  of  steel  used  and  the 
dimensional  details  of  the  heat-sensitive  elements 
or  heat  meters,  which  were  similar  to  those  de¬ 
scribed  in  reference  12.  A  sketch  showing  the 
location  of  the  heat  meters  (designated  as  stations 
1,  2,  and  3)  and  including  pertinent  dimensions 
applying  to  both  models  is  presented  in  figure  1  (a). 
For  all  the  configurations  tested,  end  plates  simi¬ 
lar  to  that  illustrated  for  the  cylinder  at  40°  yaw 
in  figure  1(a)  were  used  on  the  models.  The  dis¬ 
tances  that  the  end  plates  extended  ahead  of  the 
models  for  all  yaw  angles  are  given  in  the  figure. 
The  width  of  the  end  plates  was  held  constant  for 
all  yaw  angles. 

The  first  model  was  constructed  of  type  304 
stainless  steel,  and  water  was  used  as  the  coolant. 
This  model  had  to  be  discarded  after  several  tests 
at  40°  yaw  angle  had  been  completed  because  of 
erratic  data  obtained  as  a  result  of  the  corrosive 
action  of  the  water  on  the  electrical  leads  and  heat 
meters.  Subsequently,  a  second  model  was  con¬ 
structed  in  which  an  improved  heat-meter  design 
was  used.  The  coolant  used  in  this  model  was 
Varsol  (hydrocarbon  solvent,  boiling  point  150°F 
to  200°  F),  which  is  both  noncorrosive  and  noncon¬ 
ducting  and  may  be  cooled  to  a  lower  temperature 
than  water.  The  second  model  was  used  for  the 
remainder  of  the  tests  at  yaw  angles  of  0°,  10°,  20°, 
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Figure  1. — Schematic  diagram  of  apparatus.  All 
dimensions  are  in  inches. 


and  60°.  The  heat  meters  used  in  the  second 
model  will  be  described  in  some  detail. 

A  sectional  sketch  of  the  heat  meter  is  shown 
in  figure  1(b).  The  principle  of  operation  of  the 
heat  meter  is  described  in  detail  in  reference  12. 
Briefly,  the  two  interfaces  at  the  stainless-steel 
and  constantan  joints  are  the  hot  and  cold  junc¬ 
tions  of  a  thermocouple  circuit.  The  thermoelec¬ 
tric  output  is  then  proportional  to  the  quantity 


(b)  Enlarged  cross-sectional  sketch  of  heat  meter  used  in 
second  model. 


of  heat  flowing  through  the  constantan  disk  when 
steady  temperatures  are  attained.  A  guard  ring, 
as  shown  in  figure  1(b),  was  used  to  reduce  heat- 
conduction  errors.  Type  416  stainless  steel  was 
used  for  the  second  model  and  heat-meter  parts 
because  the  thermal  conductivity  of  this  steel  is 
approximately  equal  to  that  of  constantan  with 
the  result  that  the  meters  and  model  wall  are  more 
nearly  homogeneous.  Careful  measurements  were 
made  of  each  of  the  components  so  that  after  the 
units  were  soldered  together  the  thickness  of  the 
flanges  and  copper-constantan  joints  would  be 
known.  For  these  heat  meters  tlj,e  total  thickness 
of  the  joints  is  less  than  0.6005  inch.  In  reference 
12  it  is  shown  that  the  correction  for  heat  conduc¬ 
tion  in  the  flanges  depends  critically  on  the  flange 
thickness,  which,  for  these  meters,  was  0.001  to 
0.003  inch.  The  heat  meters  used  in  the  first  model 
were  similar  to  those  just  described  except  for  the 
flange  thickness,  which  was  about  0.008  inch,  and 
the  thermal  conductivity  of  the  304  stainless  steel 
(approximately  three-fourths  of  the  value  for 
constantan). 

Although  extreme  care  was  used  in  the  construc¬ 
tion  of  the  heat  meters  it  was  found  that  the  volt¬ 
age  outputs  of  all  three  meters  were  not  the  same 
under  similar  heating  conditions.  A  calibrator,  a 
sketch  of  which  is  shown  in  figure  1(c),  was  then 
made.  It  consisted  of  a  copper-constantan  heat 
meter  with  two  guard  rings  and  with  one  end  con¬ 
toured  to  fit  the  model.  The  second  ring  of  the 
calibrator  was  also  a  heat  meter  so  that  the  heat 
flux  through  it  could  be  measured.  The  outside 
of  the  calibrator  was  insulated.  Copper  was  used 
to  insure  a  uniform  temperature  distribution  in 
the  unit.  The  upper  end  of  the  calibrator  was 
heated  electrically.  A  0.0005-inch  plastic  film 
was  used  between  the  calibrator  and  the  model 
during  calibration  to  insulate  it  electrically  from 
the  model.  The  calibrator  was  centered  over 
a  heat  meter  and  a  constant  coolant  flow  rate 
through  the  model  was  established.  After  steady 
temperatures  were  obtained,  the  heat  flux  through 
the  calibrator  could  be  calculated  from  the  volt¬ 
age  output  and  the  known  conductivity  and  thick¬ 
ness  of  the  constantan  in  the  calibrator.  Since 
the  diameter  of  the  center  portion  of  the  calibrator 
was  the  same  as  the  heat-meter  diameter  (0.125 
inch)  and  since  the  heat  flux  through  the  second 
ring  of  the  calibrator  was  in  all  cases  approxi¬ 
mately  equal  to  the  heat  flux  through  the  center, 


Figure  1. — Continued. 
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(c)  Enlarged  cross  section  of  heat-meter  calibrator. 
Figure  1. — Concluded. 


it  was  assumed  that  this  latter  value  was  the  cor¬ 
rect  heat  flow  into  the  meter.  (Note  that  when 
the  calibrator  is  properly  centered  over  the  heat 
meter,  there  is  no  heat  flow  into  the  outer  surface 
of  the  heat-meter  flange.)  The  ratio  of  this  cal¬ 
culated  heat-flow  rate  through  the  center  of  the 
calibrator  to  the  indicated  value  from  the  heat 
meter  in  the  model  was  then  used  as  a  correction 
factor  to  apply  to  all  test  data.  The  factors  for 
the  heat  meters  in  the  second  model  were  1.22, 
0.95,  and  1.01  at  stations  1,  2,  and  3,  respectively. 

Chromel-alumel  theijmocouples  were  used  to 
measure  the  inside  surface  temperatures.  The 
leads  were  spotwelded  side  by  side  at  the  inside 
surface  of  the  heat  meters. 

The  electrical  output  of  the  heat  meters  was 
amplified  by  a  direct-current  amplifier  and  fed 
into  a  self-balancing  potentiometer.  The  output 
of  the  thermocouples  was  fed  directly  into  a  self¬ 
balancing  potentiometer.  A  similar  instrument 
was  used  to  record  stream  stagnation  temperature. 
The  stagnation  pressure  was  measured  with  a 
Bourdon  pressure  gage.  The  pressure  distribu¬ 
tions  on  the  models  were  measured  on  mercury 
manometers  whenever  possible.  Gages  were  used 
for  pressures  exceeding  60  pounds  per  square  inch 
gage.  The  locations  of  the  pressure  orifices  are 


shown  in  figure  1(a). 

The  model  was  rotated  180°  during  the  tests 
by  means  of  an  actuator  and  pulleys  so  that  a 
complete  cliordwise  survey  of  the  heat  transfer 
and  pressure  distribution  around  the  cylinder 
could  be  made  during  a  single  test.  The  surveys 
or  distributions  were  thus  obtained  in  planes 
normal  to  the  cylinder  axis.  The  temperature  of 
the  coolant  was  adjusted  by  pumping  it  through 
either  a  cooler  (alcohol  and  dry-ice  heat  exchanger) 
or  a  heater  (steam  heat  exchanger)  or  any  com¬ 
bination  of  each.  The  temperature  of  the  Varsol 
could  be  varied  from  —50°  F  to  150°  F  and  the 
water  temperature  from  60°  F  to  200°  F. 

TEST  METHODS  AND  DATA  REDUCTION 

The  tests  were  conducted  in  the  following  man¬ 
ner.  The  tunnel  and  piping  system  were  first 
preheated  to  the  desired  stagnation  temperature. 
The  model  was  polished;  then  the  temperature  of 
the  coolant  was  adjusted  and  steady  airflow  was 
established. 

Two  series  of  runs  were  made  at  each  yaw  angle. 
The  first  series  was  made  at  approximately  con¬ 
stant  tunnel  conditions  and  the  heat  rates  were 
changed  by  changing  the  model  temperature 
(varying  the  coolant  temperature).  The  results 
from  a  typical  set  of  runs  at  three  coolant  tem¬ 
peratures  are  shown  in  figure  2  where  a  heating- 

rate  parameter  is  plotted  against  a  wall 

T  —  T 

temperature  parameter  — -•  As  pointed  out 


Figure  2. — Variation  of  heat-transfer  parameter  with 
wall-temperature  parameter  at  station  2  for  typical 
runs.  A  =  60°;  250°  F;  Rm  «  1.4X  10*. 
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in  reference  12,  the  slope  of  these  lines  is  pro¬ 
portional  to  the  heat-transfer  coefficient  and  the 
recovery  temperature  may  be  obtained  from  the 
intercept  on  the  abscissa.  Sufficient  runs  of  this 
type  were  made  so  that  the  variation  of  TJT0 
with  Reynolds  number  was  obtained.  Then  the 
second  series  of  runs  was  made  to  obtain  the 
variation  of  hD/km  with  Reynolds  number.  This 
was  done  by  varying  the  stagnation  pressure 
which  ranged  from  230  to  500  pounds  per  square 
inch  gage. 

All  data  presented  in  the  report  have  been 
corrected  by  using  the  calibration  constants 
previously  discussed.  A  flange-conduction  cor¬ 
rection,  which  is  described  in  detail  in  the  appendix 
of  reference  12,  has  also  been  applied  to  all  data. 
This  correction  was  taken  to  be  a  constant  (0.7) 
for  all  the  meters  and  is  believed  to  be  correct 
to  within  ±10  percent  according  to  calculations 
made  by  the  method  of  reference  12,  which  uses 


measured  heat  rates  and  calculated  plug  tempera¬ 
tures.  The  plug  temperature  Tp  is  the  outside 
surface  temperature  calculated  from  the  measured 
inside  surface  temperature  and  the  assumption 
of  one-dimensional  heat  flow  through  the  plug. 

Radiation  has  been  neglected  since  for  all  cases 
the  radiation  is  very  small  compared  with  the 
aerodynamic  heat  transfer. 

RESULTS  AND  DISCUSSION 

PRESSURE  DISTRIBUTIONS 

Pressure-distribution  data  are  presented  because 
the}'  are  required  in  the  theoretical  heat-transfer 
calculations  and  also  to  provide  a  comparison  with 
similar  data  from  other  investigations  (for  ex¬ 
ample,  refs.  13  and  14)  where  test  conditions  and 
model  end  supports  or  configurations  are  generally 
different.  The  present  data  are  plotted  in  figure  3 
as  the  ratio  of  pi/ps  against  0*  for  A=0°,  20°, 
40°,  and  60°.  Also  shown  in  figure  3  are  curves 


Figure  3. — Pressure  distributions  on  circular  cylinder. 
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calculated  from  modified  Newtonian  theory  which 
predicts  the  pressure  distribution 

£■=/ 1  cos2  0*+ —  (1) 

ps  \  pJ  Ps 

where  pm/ps  depends  only  on  the  normal  Mach 
number  component.  The  experimental  values  of 
Pi/Ps,  however,  are  seen  to  be  essentially  inde¬ 
pendent  of  the  normal  Mach  number  or  yaw 
angle  and  would  all  fall  approximately  on  one 
curve  in  agreement  with  the  data  of  reference  13 
as  shown  in  reference  4. 

The  pressure  data  from  reference  7  have  been 
converted  to  pifps  and  are  plotted  in  figure  3. 
At  comparable  yaw  angles  and  normal  Mach 
numbers  the  data  are  in  good  agreement  except 
at  A— 60°  and  MN>m—2.07.  Thus,  even  though 
the  normal  Mach  numbers  are  about  the  same 
(2.07  in  the  present  data  and  1.98  in  ref.  7)  the 
pressures  on  the  yawed  cylinder  in  the  present 
investigation  are  less  than  those  on  the  unyawed 
cylinder  of  reference  7  up  to  about  0*=1OO°. 
The  trend  appears  to  be  the  same  as  that  shown 
in  the  data  of  reference  4  where  the  pressures 
increase  as  the  normal  Mach  number  decreases. 
This  condition  indicates  that  at  the  larger  yaw 
angles  the  pressure  distribution  on  a  yawed 
cylinder  is  not  solely  a  function  of  the  normal 
Mach  number. 

Spanwise  (parallel  to  cylinder  axis)  variations 
in  pressure  were  small  and  apparently  were 
caused  by  variations  in  tunnel  flow  which  are 
shown  in  the  tunnel  calibration  data  of  reference  7. 

The  effect  of  Reynolds  number  on  the  pressure 
distribution  is  also  small  except  on  the  back  of 
the  cylinder  in  the  separated  flow  region  where  the 
pressure  ratio  decreases  with  increasing 

HEAT  TRANSFER  AND  RECOVERY*  TEMPERATURES 

Comparison  of  data  and  theory  at  stagnation 
line, — The  values  of  the  Nusselt  number  from 
station  2  at  the  stagnation  line  of  the  cylinder 
are  plotted  against  stream  Reynolds  number  in 
figure  4  for  A=0°,  10°,  20°,  40°,  and  60°.  The 
lines  shown  in  the  figure  are  the  corresponding 
theoretical  values  calculated  for  laminar  and 
turbulent  boundary  layers. 

The  equation  used  to  calculate  the  values 
shown  in  figure  4  for  heat  transfer  in  a  turbulent 
boundary  layer  at  the  stagnation  line  of  the 


Figure  4. — Nusselt  number  at  the  stagnation  line,  span- 
wise  station  2.  Tailed  symbols  indicate  data  from 
first  model. 


cylinder  is  written 

w_(u^ykNr,,,  (^)A<sin  a)S 

[i*cosA  ® 


where  a  and  n  are  constants  in  the  Blasius  skin- 
friction  law 


For  the  present  calculations  these  values  were 
taken  as  a=0.0228  and  n— 4.  Equation  (2) 
predicts  zero  heat  transfer  at  A— 0°  and  a  peak 
in  stagnation-line  heat  transfer  at  A  ~30°  followed 
by  a  decrease  in  heat  transfer  of  about  40  percent 
from  A=40°  to  60°.  The  predicted  zero  heat 
transfer  at  the  stagnation  line  for  zero  yaw  angle 
is  not  physically  realistic  since  presumably  a 
laminar  boundary  layer  would  always  exist  for 
these  conditions  with  no  possible  mechanism 
present  to  cause  a  turbulent  boundary  layer. 
The  complete  derivation  of  equation  (2)  is  given 
in  appendix  A.  Briefly,  the  integral  momentum 
and  energy  equations  for  the  turbulent  boundary 
layer  on  a  yawed  infinite  cylinder  are  solved  by 
using  simplifying  assumptions  for  skin  friction, 
velocity  and  thermal  profiles,  secondary  flow, 
and  Reynolds  analogy.  Probably  the  most 
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restrictive  assumptions  in  this  analysis  are:  (1) 
The  flat-plat  skin-friction  law  in  the  form  of 
equation  (3)  applies  at  any  point  on  the  cylinder 
as  determined  by  local  boundary-layer  thickness 
and  local  resultant  flow  direction ;  (2)  zero  second¬ 
ary  flow  exists  over  the  entire  cylinder  (correct 
at  A=0  and  at  stagnation  line  only  for  A>0); 
(3)  and  the  flat-plate  Reynolds  analogy  applies 
between  the  heat  transfer  and  resultant  shear. 

The  equation  used  to  calculate  the  theoretical 
laminar  heat-transfer  values' shown  in  figure  4  is 

•  ^=0.77z(§f— )V2[^—(— — ITT4  (4) 

LTsPAp*  /J 

which  is  derived  in  appendix  B  and  applies  only 
for  iVp^O.7  and  Mivr,co>  1.5.  Equation  (4) 
is  based  on  the  exact  solutions  of  reference  10 
for  the  laminar  compressible  boundary  layer 
on  yawed  infinite  cylinders.  The  principal  as¬ 
sumptions  used  in  reference  10  are  constant 
specific  heats,  constant  Prandtl  number,  and 
the  perfect  gas  law.  The  variation  of  heat- 
transfer  coefficient  with  yaw  angle  as  predicted 
by  equation  (4)  can  be  closely  approximated  by 
(cos  A)1,1  for  large  stream  Mach  numbers.  (See 
appendix  B.) 

At  a  yaw  angle  of  0°  the  data  are  in  good  agree¬ 
ment  with  laminar  theory  except  at  the  lowest 
Reynolds  number  where  the  data  are  about  8 
percent  higher  than  the  theoretical  values.  At 
a  yaw  angle  of  10°  the  data  are  in  agreement  with 
values ‘obtained  by  laminar  theory  at  the  lowest 
Reynolds  number  but,  as  the  Reynolds  number 
is  increased,  the  data  become  much  higher  than 
values  indicated  by  laminar  theory  and  tend  to 
approach  the  line  calculated  by  turbulent  theory. 
The  data  at  a  yaw  angle  of  20°  appear  to  vary 
as  Rn  to  a  constant  power  (straight-line  variation) 
and  are  from  about  75  to  130  percent  higher  than 
laminar-theory  values  and  about  20  percent  less 
than  turbulent-theory  values.  At  yaw  angles 
of  40°  and  60°  the  data  are  in  good  agreement 
with  predictions  from  turbulent  theory. 

A  possible  explanation  for  the  behavior  of 
the  data  at  yaw  angles  of  10°  and  20°  is  as  follows : 
Since  the  data  .at  0°  yaw  are  in  agreement  with 
laminar-theory  calculations  and  the  data  at  40° 
and  60°  yaw  are  in  agreement  with  the  turbulent- 
theory  calculations,  it  may  be  tentatively  assumed 
that  the  theories  are  generally  valid.  Then,  since 


the  data  at  10°  and  20°  yaw  agree  with  neither 
theory,  the  boundary  layer  at  these  yaw  angles  (and 
Rm  range  of  these  tests)  is  in  some  transitional 
phase.  The  tendency  of  these  data  to  approach  a 
straight-line  variation  with  R «  (in  the  log-log  plot) 
over  such  a  large  range  of  Rm  and  A  indicates  that 
some  flow  mechanism  different  from  the  con¬ 
ventional  transitional  type  of  flow  may  be  present  . 
This  different  flow  mechanism  could  ’be  the 
vortical-type  flow  first  observed  by  the  British 
on  swept  wings  (ref.  15)  and  also  reported,  for 
example,  in  reference  16.  Another  factor  tending 
to  indicate  the  existence  of  a  flow  mechanism 
(at  yaw  angles  of  10°  and  20°)  different  from  the 
conventional  transitional  or  turbulent  boundary 
layer  is  that  the  heat  transfer  at  the  stagnation 
line  was  always  the  maximum  value  measured 
while  the  distribution  of  heat  transfer  around 
the  cylinder  decreased  monotonically  with  in¬ 
creasing  6 *  up  to  the  separation  line  of  the  chord- 
wise  flow.  This  condition  is  in  contrast  to  the 
situation  on  a  blunt  body  of  revolution  where 
conventional  transition  to  turbulent  flow  is 
generally  marked  by  a  sudden  increase  in  local 
heat  transfer  as  shown,  for  example,  in  reference 
12. 

h  D  /T  \ 

The  variation  with  yaw  angle  of  and  r 

is  shown  in  figure  5.  The  heat-transfer  data  are  a 
cross  plot  of  figure  4  at  values  of  RM  of  1.3X106, 
2.0 X 106,  and  3.5 X 106.  The  recovery-temperature 
data  are  obtained  from  plots  similar  to  figure  2  at 
Rm  ~  1 .3  X 106.  The  large  increases  in  heat  transfer 
with  increasing  yaw  angle  up  to  40°  are  clearly 
apparent  in  figure  5 .  It  may  be  noted  again  that  at 
yaw  angles  of  10°  and  20°  the  heat-transfe’r  coeffi¬ 
cients  are  below  turbulent  theory  but  at  yaw  angles 
of  40°  and  60°  good  agreement  between  data  and 
theory  is  obtained.  The  recovery  temperatures 
decrease  with  increasing  yaw  angle  as  predicted  by 
both  laminar  and  turbulent  theory  according  to 
the  equation 


(5) 


where  the  recovery  factor  r  was  taken  as  AtPt1/2 
for  laminar  boundary  layers  (see  refs.  6  and  10) 
and  NPr1/Z  for  turbulent  boundary  layers.  Al¬ 
though  the  recovery-temperature  ratios  are  less 
than  the  theoretical  values,  the  variation  with  A 
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Figure  5. — Variation  of  recovery  temperature  and  Nusselt 
number  with  yaw  angle  A  at  the  stagnation  line  of  the 
cylinder  from  station  2. 


from  A=20°  to  60°  tends  to  parallel  the  theoretical 
turbulent  variation.  The  effect  of  on  the  re¬ 
covery  temperatures  was  generally  less  than  1  per¬ 
cent  which  is  considered  to  be  within  the  range  of 
experimental  errors. 

An  indication  of  how  well  the  ideal  of  an  infinite 
cylinder  was  realized  can  be  obtained  from  figure 


6  where  the  spanwise  variations  of  are  shown 

for  A=0°,  10°,  20°,  40°,  and  60°  at  i?l  =  1.3X106, 
2X106,  and  3.5 X106.  The  data  are  cross  plots 

obtained  from  plotted  against  i?OT  for  each 

ft  QQ 


spanwise  station.  These  spanwise  stations  were 
located  2  inches  apart  along  a  generating  element 
of  the  cylinder  as  shown  in  the  sketch  of  figure  1  (a). 
Figure  6  shows  that  the  spanwise  variations- are 
within  the  expected  experimental  error  of  ±10 
percent  except  at  A=10°  and  J?„  =  1.3X106.  The 
reason  for  the  increase  in  spanwise  variations  for 
these  conditions  is  probably  that  transition  has 
just  been  initiated  at  approximately  these  values  of 
A  and  Rm  as  discussed  in  connection  with  figure  4. 


j 


(a)  i?w  =  1.3X  106. 

(b)  Ra  =2. OX  10fi. 

(c)  106. 

Figure  6. — Spanwise  variation  of  Nusselt  number  at 
stagnation  line  of  cylinder. 


Hence,  it  can  be  concluded  that,  within  the  experi¬ 
mental  accuracy,  the  heat-transfer  data  are  repre¬ 
sentative  of  conditions  on  an  infinite  cylinder. 

The  effect  of  interference  from  the  end  plates  on 
the  bow  shock  is  shown  in  figure  7  in  which  two 
schlieren  photographs  of  the  cylinder  at  a  yaw 
angle  of  60°  are  presented.  The  flow  direction  in 
these  photographs  is  from  left  to  right  and  plan- 
form  views  of  the  cylinder  and  bow  shock  are 
shown.  In  figure  7(a)  a  portion  of  the  upstream 
end  plate  is  visible  in  the  upper  left-hand  corner  of 
the  figure.  A  disturbance  which  apparently  origi¬ 
nates  at  the  leading  edge  of  the  end  plate  is  seen  to 
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intersect  the  bow  shock.  However,  the  effect  of 
this  disturbance  is  negligible  at  the  center  portion 
of  the  cylinder  since,  as  shown  in  figure  7(b),  the 
bow  shock  is  essentially  parallel  to  the  cylinder  in 
this  region. 

Circumferential  distribution  of  heat  transfer  and 
recovery  temperatures  -  and  ^  comparison  with 
theory— The  experimental  values  of  Nusselt 

number  —  at  station  2  on  the  model  are  plotted 

against  the  circumferential  angle  d*  in  figure  8  for 
A=0°,  10°,  20°,  40°,  and  60°  for  several  values 
of  rJ-  All  data  have  been  corrected  for  the  cali¬ 
bration  constant  and  heat  loss  by  conduction 
through  the  heat-meter  flange  as  discussed  pie- 
viously. 

Also  shown  in  figures  8(a)  to  8(e)  are  the  theo¬ 
retical  distributions  of  ^  for  corresponding  values 

tC  go 

of  j as  computed  for  laminar  and  turbulent 
boundary  layers  by  the  approximate  methods  dis¬ 
cussed  in  appendixes  A  and  B.  Both  the  laminar 

and  turbulent  values  of  were  calculated  from 

tC  co 

the  equation 

hD_JisD  jh 
ka~  ka  hs 

where  for  the  laminar  boundary  layer  (appendix  B) 


(a=(r 
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and  for  the  turbulent  boundary  layer  (appendix  A) 
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(a)  Flow  in  vicinity  of  (b)  Flow  in  vicinity  of 

upstream  end  plate.  station  2. 

Figure  7.— Schlieren  photographs  of  model  at  60°  yaw. 
Optical  distortion  in  schlieren  system  causes  apparent 
slight  taper  in  model.  ***  1.5X  106. 


where  the  wall  temperature  is  assumed  to  be  con¬ 
stant.  The  values  of  in  equation  (6)  were 

obtained  from  equations  (2)  or  (4).  As  shown 
in  the  appendixes,  the  only  information  needed  to 
evaluate  equations  (7)  and  (8)  is  the  pressure  dis¬ 
tribution  around  the  cylinder,  (expeflmental  pres¬ 
sure  distributions  were  used  in  all  cases  except 
at,  a=10°),  stream  Mach  number  and  stagnation 
conditions,  yaw  angle,  and  wall  temperature. 
For  the  laminar  case  the  ratios  d'J(B'w)^=i  from 
reference  10  for  2V^=1.0  are  given  in  table  I. 
The  calculations  for  all  values  of  A  were  not 
extended  beyond  6*  ~  100°  since  flow  separation 
presumably  occurred  near  this  value  of  6*  as 
indicated  by  the  approximately  constant  value 
of  pjps  for  0*>1OO°  in  the  pressure  data  of 

figure  3. 

For  a  yaw  angle  of  0°  (fig.  8(a)),  the  data  are 
seen  to  be  in  good  agreement  with  laminar  theory 
over  the  entire  forward  section  of  the  cylinder. 
The  data  (particularly  for  the  largest  Reynolds 
number)  tend  to  be  somewhat  higher  than  the 
laminar  theory  for  6*^  60°.  These  higher  values 
mav  be  caused  by  the  nonisothermal  wall  temper¬ 
atures  during  the  tests.  Both  the  laminar  and 
turbulent  calculations  from  equations  (7)  and  (8) 
are  for  an  isothermal  wall.  The  curves  for 
turbulent  theory  at  A=0°  shown  in  figure  8(a) 
were  calculated  directly  from  equations  (All)  and 
(A12)  and  are  included  for  comparative  purposes. 

The  data  at  a  yaw  angle  of  10°  (fig.  8(b))  are  in 
agreement  with  laminar  theory  at  the  lowest 
Reynolds  number;  however,  for  the  larger  value 
of  the  data  are  considerably  higher  than  the 
values  predicted  by  laminar  theory  and  up  to 
about  0*=6O°  are  still  not  so  high  as  the  values 
predicted  by  turbulent  theory.  (Both  N ewtoman 
pressure  distribution  and  experimental  ppps  from 
A=0°  Were  used  for  the  theoretical  calculations  at 
A=10°  since  experimental  pressures  were  not 
available.)  If  the  boundary  layer  at  these  larger 
values  of  were  of  the  conventional  transition 
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TABLE  I.— RATIO  OF  HEAT-TRANSFER 
COEFFICIENTS 


[From  rdf.  10;  iV>r=1.0;  linear  viscosity  relation;  zero 
transpiration  cooling] 
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0 
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1.  0574 

1.  0 

1.  0000 

1.  0000 

1.  0000 

.  5 

.  9319 

.  9118 

.  9034 

.  2 

.  8553 

.  8103 

.  7904 

0 

.  7509 

.  6510 
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type,  an  increase  in  heat  transfer  denoting  transi¬ 
tion  to  a  completely  turbulent  boundary  layer 
might  be  expected  at  some  location  downstream 
of  the  stagnation  line.  Since  this  increase  did  not 
occur  in  the  present  tests,  it  may  be  inferred  that 
a  different  flow  mechanism  is  present  such  as  the 
stream  wise  vortical  flow  observed,  for  example,  in 
references  15  and  16.  A  flow  of  this  nature  might 
be  expected  to  cause  heat-transfer  rates  inter¬ 
mediate  between  the  rates  for  a  laminar  and 
completely  turbulent  boundary  layer  and  also 
would  be  expected  to  extend  forward  to  the 
vicinity  of  the  stagnation  line  at  sufficiently  large 
Reynolds  numbers.  (See  ref.  15.)  Transition 
to  completely  turbulent  flow  has  apparently  not 
occurred  (at  least  for  75°^>0*>O)  at  A=10°  and 


(a)  A  =  0°. 

Figure  8. — Variation  of  Nusselt  number  around 
circumference  of  cylinder  at  station  2. 

3.8X106,  not  only  because  the  experimental 
values  are  below  theoretical  turbulent  values  but 
also  because  no  peak  in  heat-transfer  rates  was 
observed  as  is  predicted  b}~  the  theory  for  small 
values  of  A.  A  peak  in  the  heat-transfer  dis¬ 
tribution  is  to  be  expected  in  a  turbulent  boundary 
layer  on  cylinders  at  small  yaw  angles  by  analogy 
with  the  corresponding  theoretical  predictions  for 
an  unyawed  cylinder  or  a  body  of  revolution  as 
shown,  for  example,  in  reference  12. 

At  A=20°  (fig.  8(c)),  the  experimental  data  are 
well  above  the  theoretical  laminar  values  and 
below  the  theoretical  turbulent  values  except  for 
0*>6O°.  Here  again  it  may  be  speculated  that  a 
flow  mechanism  different  from  the  conventional 
transitional  type  may  exist  since  no  increase  in 
heat  transfer  occurred  with  increasing  0*  around 
the  cylinder.  The  fact  that  the  data  for  both 
A=10°  and  20°  are  higher  than  theoretical 
turbulent  values  for  large  values  of  0*  is  probably 
an  indication  of  errors  in  the  theory  rather  than 
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□  2.33  -  Experimental  data 
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-Modified  Newtonian  theory 


-  Experimental  ^with  -^r  from  A  =  O^j 

-  Experimental  ^with-^-  from  A  =  0°j  '“theory' 


Turbulent 

theory 


(b)  A=  10°. 

Figure  8. — Continued. 


the  possible  development  of  a  turbulent  boundary 
layer  on  the  cylinder.  That  is,  the  errors  caused 
by  the  simplifying  assumptions,  such  as  zero 
secondary  flow,  isothermal  wall,  and  flat-plate 
skin  friction  are  expected  to  increase  with  increas¬ 
ing  distance  from  the  stagnation  line.  (See 
appendix  A.) 

In  contrast  to  the  results  at  A— 10°  and  20°, 
the  data  at  A=40°  (fig.  8(d))  are  in  fairly  good 
agreement  or  are  generally  somewhat  higher  than 
the  theoretical  turbulent  distributions  over  the 
entire  forward  region  of  attached  flow.  At  A=60° 


0  20  40  60  80  100  120  140  160  180 

0*  deg 

(c)  A  — 20°. 

Figure  8. — Continued. 

(fig.  8(e))  the  theory  underestimates  the  heat- 
transfer  coefficient  for  large  values  of  0*  by  as 
much  as  50  percent  whereas  at  A— 40°  the  agree¬ 
ment  between  data  and  theor\~  is  good  except  at 
i?co  =  1.49X106  and  2.75X10®  with  the  second 
model.  On  the  basis  of  this  comparison,  the 
boundary  layer  is  probably  completely  turbulent 
(including  the  flow  at  the  stagnation  line)  for  yaw 
angles  of  about  40°  or  greater.  The  larger  dis¬ 
crepancy  between  theory  and  data  from  the  second 
model  at  A=40°  is  believed  to  be  caused  by  differ¬ 
ences  in  heat-meter  construction  as  discussed  in 
the  section  entitled  “Apparatus.”  The  data  from 
the  second  model  are  believed  to  be  more  nearly 
correct. 

Downstream  from  the  separation  line,  which 
from  figure  3  is  in  the  vicinity  of  0*  =  11O°,  the 
heat-transfer  coefficients  for  all  values  of  A  are 
approximately  constant.  The  values  of  the  heat- 
transfer  coefficients  in  this  separated-flow  region 
are  from  about  10  to  15  percent  of  the  stagnation- 
line  value  except  at  A=0°  and  10°  with  Ra 
<2X106. 

Experimental  and  theoretical  values  of  the  ratio 
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of  local  recovery  temperature  to  stagnation  tem¬ 
perature  Tej T0  are  shown  in  figure  9.  The  experi¬ 
mental  data  were  obtained  from  plots  similar  to 
figure  2  and  the  theoretical  values  were  calculated 
from  equation  (5)  after  replacing  Ts/T0  with 
TijT0.  Experimental  values  of  px lps  were  used  to 
evaluate  TJT0  except  for  A=10°  where  the  pres¬ 
sure  ratio  distribution  for  A— 0°  was  used.  Over 
the  front  half  of  the  cylinder  the  recovery- 
temperature  data  are  always  below  those  for  the 
turbulent  theory  (figs.  9(a)  to  9(e)),  and  at 
A=0°,  10°,  and  40°,  the  data  are  also  somewhat 
below  the  laminar  theory  but  at  A=20°  the  data 
tend  to  approach  the  turbulent  theory  for  0*>6O°. 
At  A=60°  (fig.  9(e)),  the  data  are  between  the 
curves  for  the  laminar  and  turbulent  theory.  In 
general,  then,  no  definite  indipalions^as  to  which 
type  of  flow  is  present  can  be  obtained  from  the 
recovery-temperature  data;  however,  the  small 


(a)  A  =  0°. 

(b)  A=  10°. 

(c)  A  =  20°. 

Figure  9. — Variation  of  recovery  temperature  around 
circumference  of  cylinder  at  station  2. 
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(d)  A  =  40°. 

(e)  A  =  60°. 

Figure  9. — Concluded. 


differences  involved  are  probably  within  the 
experimental  error. 

At  A  ^=40°  (fig.  9(d))  the  values  of  TJT0  for  the 
first  model  are  smaller  for  0*>3O°  than  the  values 
from  the  second  model.  This  discrepancy  is 
again  presumably  caused  by  differences  in  heat- 
meter  construction. 

Effect  of  chordwise  temperature  gradient. — The 

differences  in  the  heat  transfer  and  recovery 
temperature  distributions  between  the  first  and 
second  models  at  A =40°  as  shown  in  figures 
8(d)  and  9(d)  may  have  been  caused  partly  by  the 
different  chordwise  wall-temperature  distributions 
for  the  two  models.  (Chordwise  is  here  used  in 
the  sense  of  describing  the  circumferential  dis¬ 
tribution  in  a  plane  normal  to  the  cylinder  axis.) 
These  wall  temperature  distributions  are  different 
because  of  the  different  coolants  used,  water  in 
the  first  model  and  Varsol  in  the  second  model. 
An  attempt  was  made  to  evaluate  this  effect  of 
wall  temperature  gradients  from  an  additional  set 
of  tests  with  the  second  model  using  water  as  a 
coolant  instead  of  Varsol. 


In  figure  10  the  ratio  TpjT0  is  plotted  against 
8 *  for  typical  runs  with  the  second  model  for  which 
water  and  Varsol  were  used  as  coolants.  Since 
water  is  a  more  efficient  coolant,  the  temperature 
distribution  for  the  water-cooled  model  is  more 
nearly  isothermal  than  that  for  the  Varsol  cooled 
model,  and  hence  the  chordwise  temperature 
gradients  are  smaller  for  the  water-cooled  model. 

Also  shown  in  the  upper  part  of  figure  10  are 
the  ratios  of  equilibrium  to  stagnation  tempera¬ 
tures  that  were  obtained  by  using  the  two  coolants. 
A  qualitative  analysis  of  the  effect  of  upstream 
temperature  distributions  such  as  those  of  figure 
10  on  local  equilibrium  temperature  (wall  tem¬ 
perature  required  for  zero  heat  transfer)  indicates 
that,  as  the  temperature  distribution  becomes 
more  nearly  isothermal,  the  local  equilibrium 
temperature  would  decrease.  This  result  is  in 
agreement  with  the  present  data  since  the  water- 
cooled  model  was  always  more  nearly  isothermal 
than  the  Varsol  cooled  model. 

It  can  also  be  shown  that  the  local  heat  rate 
would  be  smaller  when  the  upstream  temperature 
distribution  is  more  nearly  isothermal  if  the  local 
temperature  potential  (difference  between  recov¬ 
ery  and  wall  temperature)  is  held  constant. 

The  effect  of  wall  temperature  gradients  on  the 
Nusselt  number  or  heat-transfer  coefficient  is 
determined  by  the  relative  changes  in  both  the 
heat  rates  and  temperature  potential.  Nusselt 
number  distributions  from  the  tests  on  the  second 


Figure  10. — Effect  of  coolant  on  temperature  distribu¬ 
tions  for  A  =  40°  and  =  2.75X  106.  (Data  from 
station  1  on  the  second  model.) 
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model  with  water  and  Varsol  used  as  coolants  are 
shown  in  figure  11.  Comparison  of  these  data 
indicates  that  the  heat-transfer  coefficient  de¬ 
creases  more  with  increasing  0*  on  the  water-cooled 
model. 

The  data  of  figure  10  were  used  to  calculate  the 
effect  of  wall  temperature  gradients  in  a  turbulent 
boundary  layer  by  the  method  of  reference  17. 
This  method  is  based  on  an  approximate  solution 
of  the  integral  energy  equation  for  the  incompres¬ 
sible  turbulent  boundary  layer  on  a  flat  plate  with 
a  stepwise  temperature  distribution  at  the  surface 
of  the  plate.  The  method  is  extended  approxi¬ 
mately  to  the  case  of  a  pressure  gradient  by  using 
the  difference  Tp~  Te  rather  than  Tv  as  the  depend¬ 
ent  variable.  The  results  of  the  calculation  for 
the  temperature  distributions  of  figure  10  indi¬ 
cated  that  at  0*=9O°  the  heat-transfer  coefficient 
on  the  Varsol  cooled  model  would  be  about  11 
percent  larger  than  the  value  on  the  water-cooled 
model.  This  value  compares  with  a  value  of  20 
percent  from  the  data  of  figure  11. 

This  discussion  indicates  that  the  differences 
in  the  data  as  shown  in  figure  11  may  be  partly 
caused  by  the  different  wall  temperature  gradients. 
However,  the  effects  are  within  the  experimental 
accuracy  so  that  no  definite  conclusion  is  possible 
except  that  the  larger  differences  between  the  data 
from  the  first  and  second  models  are  not  caused 
entirely  by  wall  temperature  gradients.  It  should 
be  noted  that,  although  local  heat-transfer  coeffi¬ 
cients  at  large  values  of  may  be  influenced  by 
temperatffre-gradient  effects,  the  average  or  total 
heat  transfer  to  the  cylinder  would  be  changed  at 


Figure  11. — Effect  of  coolant  on  Nusselt  number  distri¬ 
bution  at  station  1  on  second  model.  A  =  40°; 

—  2.75X  106. 


most  by  only  a  few  percent  since  most  of  the  total 
heat  is  transferred  at  small  values  of  0*. 

Average  heat  transfer  over  front  half  of  cylin¬ 
der. — The  average  Nusselt  number  over  the 
forward  portion  of  the  cylinder  from  0*=O°  to  90° 
was  obtained  by  integration  of  plots  similar  to 
figure  8.  The  results  have  been  plotted  against 
Rm  in  figure  12.  Also  plotted  in  the  figure  for 
comparison  are  data  from  reference  7  where  the 
average  heat  transfer  to  the  front  half  of  a  K-inch- 
diameter  copper  cylinder  was  measured  directly. 
At  comparable  yaw  angles  the  two  sets  of  data  are 
seen  to  be  in  good  agreement.  Since  the  small 
copper  model  of  reference  7  presumably  had 
essentially  an  isothermal  wall  temperature  distri¬ 
bution,  this  good  agreement  indicates  that  the  wall 
temperature  variations  encountered  in  the  present 
tests  had  very  little  effect  on  the  average  heat- 
transfer  coefficients.  The  reason  for  this  small 
effect  can  be  seen  from  figure  11  which  shows  that 
different  wall  temperature  gradients  such  as  those 
of  figure  10  could  not  have  any  significant  effect 
on  the  average  heat-transfer  coefficient  over 
the  front  half  of  the  cylinder.  It  therefore  seems 
reasonable  that  the  average  heat-transfer  coeffi¬ 
cients  for  an  isothermal  model  would  be  about  the 
same  as  the  present  results  within  the  experimental 
error. 

Also  shown  in  figure  12  are  the  average  heat- 
transfer  coefficients  from  the  turbulent  and  lami¬ 
nar  theories  of  the  present  report  as  discussed  in 
appendixes  A  and  B.  At  A— 0°  the  data  are  in 
good  agreement  with  laminar  theory  but  at  A=10° 
and  20°  (fig.  12(a))  the  data  increase  with  in¬ 
creasing  Keynolds  number  and  tend  to  approach 
the  turbulent  theory.  Note  that  these  average 
data  are  generally  in  better  agreement  with  the 
corresponding  theoretical  results  than  the  stag¬ 
nation-line  data  of  figure  4.  The  reason  for  this 
better  agreement  is  apparent  from  the  distribu¬ 
tions  shown  in  figures  8(b)  and  8(c).  At  A— 40° 
(fig.  12(b))  the  data  are  in  good  agreement  with 
turbulent  theory.  At  A  — 60°  the  data  are  about 
15  percent  higher  than  the  present  theoretical 
values  but  are  in  good  agreement  with  the  theoret¬ 
ical  values  of  reference  7.  This  latter  theory  is 
based  on  an  average  flat-plate  correlation  formula 
using  a  reference  temperature  and  the  cross-flow 
velocity  for  a  yawed  cylinder.  The  value  of  an 
empirical  constant  was  obtained  from  the  data  of 
reference  7.  The  agreement  between  the  present 
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(a)  A=0°,  10°,  and  20°. 

Figure  12.  Average  Nusselt  number  on  front  half  of 
circular  cylinder. 


/Poo 

(b)  A  =  40°  and  60°. 
Figure  12. — Concluded. 


data  and  the  data  of  reference  7  indicates  that 
this  same  constant  also  applies  to  the  present 
data.  The  present  theory  for  a  turbulent  bound¬ 
ary  layer  (appendix  A)  is  apparently  in  consider¬ 
able  error  at  large  values  of  d*  when  A  is  also  large, 
possibly  because  of  the  simplifying  assumption  of 
zero  secondary  flow  which,  for  large  values  of  6*, 
would  tend  to  increase  with  increasing  A.  This 
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error  is  also  apparent  in  figure  8(e)  which  shows 
that  the  present  theory  underestimates  the  local 
heat  transfer  by  50  percent  at  0*=9O°  and  A=60°. 

CONCLUDING  REMARKS 

Local  heat  transfer,  equilibrium  temperatures, 
and  wall  static  pressures  have  been  measured  on 
a  circular  cylinder  at  yaw  angles  of  0°,  10°,  20°, 
40°,  and  60°.  The  tests  were  made  at  a  stream 
Mach  number  of  4.15  and  stream  Reynolds  num¬ 
ber  range  of  1  X106  to  4  X  106  based  on  cylinder 
diameter. 

The  pressure  distributions  up  to  the  separation 
line  of  the  chordwise  flow  •werevessentially  inde¬ 
pendent  of  normal  Mach  'number  and  yaw  angle 
A.  *  6 

At  the  stagnation  line  of  the  cylinder  the  heat- 
transfer  coefficients  increased  from  100  to  180 
percent,  depending  on  the  Reynolds  number,  as 
the  yaw  angle  was  increased  from  0°  to  40°.  In¬ 
creasing  the  yaw  angle  from  40°  to  60°  resulted 
in  a  40-percent  decrease  in  heat-transfer  coeffi¬ 
cients.  At  zero  yaw  the  data  are  in  agreement 
with  theory  for  laminar  boundary  layers  and  at 
40°  and  60°  yaw  the  data  agree  with  a  theory  for 
turbulent  boundary  layers.  At  yaw  angles  of  10° 
and  20°  the  data  were  generally  somewhere  be¬ 
tween  the  theoretical  laminar  and  turbulent  values. 

The  chordwise  distribution  of  heat-transfer  co¬ 
efficients  up  to  the  separation  line  was  |n  good 
agreement  with  laminar  theory  at  zero  yaw  angle. 
At  yaw  angles  of  10°  and  20°  the  chordwise  dis¬ 
tributions  were  generally  above  laminar  theory 
and  below  turbulent  theory.  At  any  given  yaw 
angle  the  heat  transfer  at  the  stagnation  line  was 
the  peak  value  and  the  chordwise  distribution  al¬ 
ways  decreased  monotonically  from  this  peak 
value.  The  values  of  the  heat-transfer  coefficients 
as  well  as  the  chordwise  distributions  at  yaw  angles 
of  10°  and  20°  indicated  that  a  boundary-layer 
flow  mechanism  different  from  the  conventional 
transitional  or  turbulent  flow  may  have  been 
present.  At  yaw-  angles  of  40°  and  60°  the  chord- 
wise  distributions  were  in  good  agreement  with 
turbulent  theory  except  for  large  values  of  0* 
(angular  distance  from  the  stagnation  line)  at 
A=60°.  On  the  downstream  side  of  the  cylinder 
in  the  separated-flow  region  where  the  theories  do 
not  apply,  the  heat-transfer  coefficients  were  ap- 
pioximately  constant  at  about  10  to  15  percent  of 
the  stagnation-line  value. 
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The  average  heat-transfer  coefficients  over  the 
front  half  of  the  cylinder  are  in  agreement  with 
previous  data  for  a  comparable  Reynolds  number 
range. 

The  theoretical  heat-transfer  distributions  for 
both  laminar  and  turbulent  boundary  layers  are 
calculated  directly  from  simple  quadrature  for¬ 
mulas  derived  in  the  appendixes.  Simple  formulas 
for  calculating  the  stagnation-line  heat  transfer 
with  both  laminar  and  tiirbulent  boundary  layers 


are  also  presented.  The  theoretical  variation  of 
stagnation-line  heat-transfer  coefficient  with  yaw 
angle  A  is  approximately  (cos  A)1*1  for  laminar 
boundary  layers  but  for  turbulent  boundary  layers 
the  heat-transfer  coefficient  increases  at  first, 
reaches  a  maximum  at  A«30°,  and  then  decreases 
about  40  percent  from  A=40°  to  60°. 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Langley  Air  Force  Base,  Va.,  December  2 ,  1958. 


APPENDIX  A 


TURBULENT  BOUNDARY-LAYER  THEORY  ON  YAWED  INFINITE  CYLINDERS 


The  starting  point  of  the  present  analysis  is  the 
integral  equations  of, momentum  and  energy  in 
the  compressible  boundary  layer  on  a  yawed  in¬ 
finite  cylinder.  These  equations  may  be  obtained 
directly  from  appendix  B  of  reference  10  by  re¬ 
placing  the  expressions  for  laminar  shear  and  heat 
transfer  with  general  expressions  which  include 
the  effects  of  turbulent  shear  and  heat  transfer. 
In  terms  of  a  modified  Stewartson  coordinate 
system  and  for  zero  normal  velocity  at  the  wall, 
these  equations  become 
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The* transformed  coordinates  and  velocities  are 
related  to  the  corresponding  physical  quantities 
by  the  following  equations: 


The  shear  and  heat  transfer  in  the  transformed 
system  are  related  to  the  corresponding  physical 
quantities  by  the  following  definitions: 


T 


PoPo  To 
PrPr  Tx 


(A5) 


where  the  subscript  r  denotes  quantities  evaluated 
at  some  reference  point  within  the  boundary 
layer.  Essentially,  the  same  transformation  as 
equations  (A4)  has  been  used  in  reference  18  where 
good  correlation  with  experimental  skin-friction 
data  on  insulated  plates  was  obtained  by  taking 
the  reference  point  in  the  local  free  stream.  In 
reference  19  the  transformation  was  used  to  cal¬ 
culate  heat  transfer  in  turbulent  boundary  layers 
on  blunt  bodies,  and  better  agreement  with  ex¬ 
perimental  data  was  generally  obtained  by  evaluat¬ 
ing  the  reference  quantities  so  that  the  flat-plate 
skin-friction  relation  of  reference  20  was  satisfied 
when  the  method  was  applied  to  this  case. 

In  order  to  get  a  practical  solution  to  this 
system  of  equations,  several  simplifying  assump¬ 
tions  are  made.  The  basic  ideas  of  the  “line-of- 
flow  principle”  as  discussed  in  reference  21  are 
adopted  in  contradiction  to  the  “independence 
principle”  of  reference  22.  The  data  of  reference 
23,  as  well  as  the  discussion  of  reference  21,  indicate 
that  the  independence  principle  does  not  apply  to 
the  growth  of  a  turbulent  boundary  layer  on  a 
yawed  flat  plate.  In  accordance  with  the  line-of- 
flow  principle,  it  is  assumed  that  the  local  resultant 
shear  is  a  function  of  the  resultant  velocity  at  the 
edge  of  the  boundary  layer  and  that  this  shear  is 
related  to  the  heat  transfer  by  a  Reynolds  analogy. 
The  effects  of  a  chordwise  pressure  gradient  are 
accounted  for  approximately  by  expressing  the 
shear  in  terms  of  local  boundary-layer  thickness. 
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The  effects  of  compressibility  are  accounted  for  by 
expressing  the  skin-friction  law  in  the  transformed 
system.  The  equation  for  the  transformed  skin 
friction  is  then  assumed  to  be  of  the  Blasius  type 
and  is  written  as 


!tx*+t  2 


=poQi2a(^r  (Ae) 


where  8  is  the  velocity  boundary-layer  thickness 
and  Q\—^Ui2jrVi  is  the  local  transformed  re¬ 
sultant  velocity  at  the- edge  of  the  boundary  layer. 
The  Colburn  modification  (ref.  *24)  of  the  Reynolds 
analogy,  also  expressed  in  the  transformed  system, 
is 

_ _ =—\  -  Tr  (A7'\ 

Cpp0Q\(Taw-~Tw)  NP2«  p0Qi2 

In  general,  a  secondary  flow  would  exist  and  its 
limiting  direction  aw  at  the  w&l\  with  respect  to 
the  cylinder  axis  is  determined  by  the  ratio  of 
chordwise  to  spanwise  shears  as  given  by 


TX 

tan  aw=— 


The  velocity  profiles  are  now  assumed  to  be 


u_(i 

7\  1/7 

ut  \  < 

V 

V  (i 

5\1/5 

F~V< 

5/ 

so  that,  from  the  definition  of  E  and  the  assump¬ 
tion  that  8X—8Y, 

E_49  ,Ainv 

5  376  ^A10^ 

Then,  from  equations  (A6),  (A8),  and  (A.10), 
equation  (A2)  may  be  written  as 

^  (tt  st\ _ 376  fv<)Ui\1/nQ\2  cos  aw  /A1-.N 

dx  a  {-&)  ft  WW*  (All) 

which,  for  the  required  boundary  condition  of 
finite  8  at  X~§,  integrates  to 

3/6  1  -\~7i 

U  i  w 

(A  12) 

From  equations  (A5),  (A6),and  (A7)  the  general 


expression  for  the  heat-transfer  coefficient  is 

1  —  _  a  ( J!o\1/U  r  Q  PrPr  IT i 

n-Taw-T-NP»*  KQj)  c»pMl  poPo  V  T0 

'  (A13) 

HEAT  TRANSFER  AT  THE  STAGNATION  LINE 

On  a  cylinder  with  a  blunt  leading  edge  U^O, 
Q\-*Vi,  and  aw~>0  as  X— >0;  hence,  from  equation 
(All),  the  expression  for  8S  is 


1±»  376  (Vx  n 

«•  ”  -wa^lln(^h 


At  the  stagnation  line,  equation  (A  13)  becomes 


—  I ~Ts(PtPt 


Tq  \PoPo 


By  substituting  equation  (Al4)  for  with 

(VOx=o=z^YUR’m  sin  A 

/ dU A  / PoPo\  T0  / duj\ 

\dX/Xs. o  \PrPr/$  Ts \dx  ) s 

(from  equations  (A4))  equation  (A  15)  may  be 
written,  after  some  rearranging,  as 


hsl  -xj  1/3/ Ur,  00  Lf-w 


(a  wry 

\  P0P0J  1 


(sin  A)  W+1  ^7 


n—1  1 

n+ 1  /  49  l  duA  i+» 

\376COS  ^u2Hx) . 


(Ate) 


In  terms  of  a  stream  Nusselt  and  Reynolds 
number,  equation  (A  16)  becomes 

n  n 

hi _  ~kt  1/3  / Ur,  co l\  (  Ps  Pr  TX  H-n 


= NPru3  1+n(aJk»L  X 

V  ^co  /  V  Pay  Po  d  T  / 

11—  1 

/  •  , \n+’/  49  n0  l  du j 

(sin  A)  I  ;r—  —  cos  A - j- 

\376ai„  dx 


l  du{\  i+» 


)r  (Ait) 


for  constant  values  of  cv  and  NPr.  The  reference 
quantities  pr  and  Tr  were  evaluated  at  the  wall 
temperature  for  comparison  with  the  experimental 
data  of  the  present  report. 
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CHORDWISE  DISTRIBUTION  OF  HEAT  TRANSFER  where 


According  to  equations  (A  13)  and  (A  14)  the 
chordwise  distribution  of  heat-transfer  coefficient 
would  be 


<1 T  (w)'' Yf.JJF,  (A18) 


If  the  flow  at  the  edge  of  the  boundary  layer  is 
isentropic,  then  the  rel a tioiy  between  temperature 
and  pressure  is 

s  \pj 


The  adiabatic  energy  equation  requires  that 
uf+vf  2 


d0 


tO -t) 


so  that  from  the  definitions  of  the  transformed 
velocities 


Hence, 


/Ts^Y2 
_Ql=Ql_1T1  To  I 
I  i  _T,J 


and  equation  (Al8)  may  be  written  as 
* 

7-1 


37-l 

hs  f^r,$TT  \Ps/ 


(  P±\~_T, 

W  To 


1- 


,n— 1 
2  n 


(!) 


1  In 


(A19) 


From  equations  (A12),  (A14),  and  (A4),  the  ratio 
of  boundary-layer  thicknesses  is 


/  5 \1/n 

W 


/TofduA 

l+n\  Ts\dx  Js  CX  VrTr.s  JT  1/rt 

Jo  1 


(£f(0 


37-1 

27  COS  Ciu 


COS  ax 


dx 


1 

1  +  n 


.(A20) 
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Va 

'GH 

The  only  quantity  as  yet  unspecified  in  equation 
(A20)  is  aw,  the  angle  of  the  “wall  streamline” 
with  respect  to  the  cylinder  ajis.  This  angle 
could  be  determined  from  the  simultaneous  solu¬ 
tion  of  equations  (Al),  (A2),  and  (A3)  by  using 
the  shear  and  heat-transfer  assumptions  of  equa¬ 
tions  (A6)  and  (A7).  Such  a  solution  is  difficult 
to  justify  in  view  of  the  lack  of  knowledge  concern¬ 
ing  turbulent  boundary  layers  on  yawed  bodies. 
Rather,  the  simplifying  assumption  of  zero  sec¬ 
ondary  flow  is  introduced  at  this  point.  That  is 


«„=«!= tan 


,  Ui 

v, 


(A21) 


which  is  exactly  correct  on  the  stagnation  line 
and  is  expected  to  be  a  good  approximation  in 
the  vicinity  of  the  stagnation  line.  The  second¬ 
ary  flow  affects  the  heat  transfer  only  through 

COS  Oi 

the  factor - -  which,  for  n— 4,  enters  equation 

(A  19)  to  the  1/5  power.  This  small  power  indi¬ 
cates  that  the  assumption  of  zero  secondary  flow 
may  not  have  much  effect  on  the  heat-transfer 
calculation  except  at  large  values  of  X  or  near 
a  separation  line.  Evaluating  the  reference  quan¬ 
tities  at  the  wall  and  assuming  that  the  wall  is 
isothermal  then  gives  the  form  of  h/hs  shown  in 
the  text  as  equation  (8)  where  8s/8  is  obtained 

from  equation  (A20)  with  C0S  ^=1  and  ^r^W-=l. 

From  equations  (A7),  (A8),  and  (A21),  the 
assumption  of  zero  secondary  flow  results  in  the 
following  analogy  between  heat  transfer  and 
spanwise  shear 


<lw  1  rr 

cPPoVi(T aw—T v)~~ N Prm  p0V 2 
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or 

T  y  Vll  _  2 

rr2  k0(Talc-T,c) 

2p°Vl 

This  relation  may  be  written  in  terms  of  physical 
shear,  heat  transfer,  and  velocities  as 


which  is  of  the  same  form  as  the  corresponding 
analogy  between  heat  transfer  and  spanwise  shear 
for  the  exact  laminar  solutions  of  reference  6.  In 
terms  of  the  chordwise  shear,  equation  (A22) 
becomes  the  flat-plate  Reynolds  analogy 

TX  PoUyl _  2  hi 

1  2  Mo  NpTm  k0 

o  Po^l 


ty  PoVil  2  hi  /aoo\  which  indicates  that  for  zero  secondary  flow  the 

1  g  Mo  ~NPrv?  To  \  ^  )  effects  of  chordwise  pressure  gradient  on  Reynolds 

2  P°Vl  .  analogy  are  neglected. 


APPENDIX  B 

HEAT  TRANSFER  IN  THE  LAMINAR  BOUNDARY  LAYER  ON  YAWED  INFINITE  CYLINDERS 


SIMPLIFIED  SOLUTION  AT  THE  STAGNATION  LINE 

The  results  of  the  exact  solutions  of  reference  10 
are  used  to  calculate  the  heat  transfer  in  a  laminar 
boundary  layer.  Thje  basic  equation  for  the  heat- 
transfer  coefficient  at  the-  stagnation  line  of  a 
yawed  infinite  cylinder  is  * 

<bi> 

where  for  NPr=0.7,  h  is  given  approximately  by 
the  correlating  equation 

£=0.5  (B2) 

\PwVwS 

A  general  expression  for  the  stream  Nusselt  num¬ 
ber  is  then 


Substitution  of  equation  (B4)  into  equation 
(B3)  then  results  in  the  equation 


(B5) 

where  it  has  been  assumed  that 


Because  of  the  small  value  of  the  exponent,  this 
assumption  is  correct  to  within  1  or  2  percent  over 
a  wide  range  of  temperatures.  Equation  (Bo) 
applies  only  for  NPr= 0.7  and 

The  effect  of  yaw  angle  on  the  heat-transfer 
coefficient  is  obtained  from  equation  (B5)  as 


jh ^ _ q  ej  I P co^R,  co  l  k w 

k  CO  Y  Moo  K 


J—  —  cos  A  (—  ( ^Y" 

Y  Pw  p co  \um  dx  J s  \picPid/ 


which,  for  constant  cv  and  NPr,  becomes 


hsD  ^-r  /  D  du{\  “|1/2 

L  V^co  dx) J 


(Psps 

\PsJk  W  ) 


(  EiEi  y/2 

VPoo/tco/ 


(B3) 


in  terms  of  the  diameter  D. 

The  stagnation-line  velocity  gradient  calculated 
from  the  modified  Newtonian  pressure  gradient  is 


which  according  to  reference  6  is  in  good  agree¬ 
ment  with  experimental  data  for  MN;m>  1.5. 
For  a  perfect  gas,  this  expression  for  the  velocity 
gradient  becomes 


/  D  du{\  _  2 

W  dx  )s~Mn>  m 


1/2 

(B4) 


which  for  a  perfect  gas  and  Sutherland's  viscosity 
temperature  law  may  be  closely  approximated  by 
(cos  A)11  for  large  free-stream  Mach  numbers. 
Equation  (B6)  shows  that  for  given  values  of 
and  A  the  ratio  hs/(hs) a=o  increases  with 
an  increase  in  temperature  level  because  of  the 
increase  in  p$/p0.  This  increase  indicates  that 
the  effect  of  yaw  as  measured  in  a  wind  tunnel 
may  be  somewhat  different  than  that  measured 
in  flight.  Thus,  for  Mm  =  4.0  and  A=60°, 
hj(hs) a=o  is  about  9  percent  larger  at  Tm  — 
390°  R  than  at  ^  =  160°  R.  It  is  of  interest  to 
compare  equation  (B6)  with  the  result  of  the 
approximate  analysis  of  reference  4.  Thus,  if  a 
linear  viscosity-temperature  relation  is  assumed 
and  the  restriction  of  large  normal  Mach  number 
(a  small  yaw  angle  being  implied)  is  imposed  so 


k  fTsyvpX1* 

(^$)a=o  \T0/  \Pt/ 
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which  is  exactly  the  same  as  equation  (11)  of  Solving  for  0  gives 
reference  4. 

CHORDWISE  VARIATION  OF  HEAT-TRANSFER  COEFFICIENT 

A  general  expression  for  the  heat-transfer  rate 
from  reference  10  is 


dX 

p=zUx  ,  dUj 
X+dX 


where  y  is  a  similarity  variable  which  can  be 
written  as 


The  ratio  of  the  local  heat  transfer  to  the  value 
at  the  stagnation  line  for  isentropic  flow  at  the 
edge  of  the  boundary  layer  is  then 


The  simplest  way  to  evaluate  0  is  to  calculate 
Ui  and  X  from  the  required  pressure-distribution 
data  by  means  of  equations  (B9)  and  then  sub¬ 
stitute  these  values  directly  into  equation  (Bll). 
(This  procedure  is  analogous  to  that  of  ref.  25.) 
The  values  of  Ui  used  in  the  present  calculations 
are  shown  plotted  against  X/D  in  figure  13  which 
also  includes  for  comparison  the  values  calculated 
from  Newtonian  pressure  distribution.  The  heat- 
transfer  coefficient  distribution  then  follows  from 
equation  (B8)  as 
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where  for  the  present  calculation  it  is  assumed 
that 


rt) 

\hs/rfPr?s  j.o  yls/ xPr=  1 .( 


where 

Hx-Hk 


t  \T o’  J 

'l~'X2\  Ul~a°\ ; 

y__  iTi 

Jo  Po  Po  M  t0 
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If  it  is  assumed  that  6  is  a  function  only  of  v,  then 
the  similar  solutions  of  reference  10  may  be  used 
where  6  has  been  calculated  for  various  values  of 
TJT0,  T0/Ts,  and  fi  with  NPr—l.  The  ratio 
T0/Ts  is  a  parameter  that  depends  on  yaw  angle 
and  Mach  number  and  0  is  the  velocity-gradient 
parameter  defined  by  the  velocity  distribution  for 
a  similar  flow  as 


Then,  since  for  NPr=1.0,  Taw=T0,  the  following 
equation  may  be  written: 

(£)  _(4k)  _r  I 

\^s/Ifpr=l.O  */Npr=1.0  L(^m)/J=1.oJ Vpr=1  0 

The  quantity  d'a  is  tabulated  in  reference  10  as 
a  function  of  TJT0,  T0/Ts,  and  0  for  NPr=  1. 
Values  of  the  ratio  6’J(6’W)^X  are  given  in  table 
I  in  the  present  report. 

An  alternative  procedure  analogous  to  the 
method  of  reference  26  is  one  in  which  an 
“equivalent”  wedge  flow  is  found  for  each  segment 
of  the  given  velocity  distribution.  The  value  of 
0  is  found  from  the  simultaneous  solution  of 
equation  (BlO)  evaluated  at  the  end  points  of 
each  segment  and  the  equation  (from  ref.  10) 


r  vnI/2 

©=|_(2— 9’u 


Ux=CX 


(BIO) 


evaluated  at  the  first  point  of  the  segment  where 
X  is  now  considered  as  the  chordwise  length 
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Figure  13. — The  variation  of  transformed  chordwise  velocity  with  transformed  chordwise  coordinate  according  to  experi- 

mental  and  Newtonian  pressure  distribution. 


from  the  leading  edge  of  the  equivalent  wedge 
flow.  A  sample  calculation  carried  out  by  this 
procedure  resulted  in  a  heat-transfer  coefficient 
about  7  percent  less  at  0*=5O°  and  A=0°  than 
the  value  obtained  by  using  the  first  procedure. 
The  method  of  reference  26  also  resulted  in 


larger  values  of  0,  namel}^  0—2.00  at  0*=51.5° 
compared  with  0=1.38  obtained  by  using  the 
method  of  reference  25.  Since  0—2  is  the 
maximum  value  (from  eq.  (Bll))  in  the  similar 
solutions,  the  method  of  reference  26  could  not 
be  used  for  0*>51.5°  in  this  particular  example. 
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